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Abstract 

In this paper, we prove a uniqueness result in the inverse problem of determining several non- 
constant coefficients of one-dimensional reaction-diffusion equations. Such reaction-diffusion 
equations include the classical model of Kolmogorov, Petrovsky and Piskunov as well as more 
sophisticated models from biology. When the reaction term contains an unknown polynomial 
part of degree N, with non-constant coefficients ^k{x), our result gives a sufficient condition 
for the uniqueness of the determination of this polynomial part. This sufficient condition only 
involves pointwise measurements of the solution u of the reaction-diffusion equation and of 
its spatial derivative du/dx at a single point xq, during a time interval (0,e). In addition to 
this uniqueness result, we give several counter-examples to uniqueness, which emphasize the 
optimality of our assumptions. Finally, in the particular cases N = 2 and = 3, we show that 
such pointwise measurements can allow an efficient numerical determination of the unknown 
polynomial reaction term. 

Keywords: reaction-diffusion • inverse problem • multi-parameter • heterogeneous media • uniqueness 

1 Introduction 

Reaction-diffusion equations arise as models in many fields of mathematical biology [20j. From 
morphogenesis [33] to population genetics [HI [T7] and spatial ecology [271 ISSl [32] > these partial 
differential equations benefit from a well-developed mathematical theory. 

In the context of spatial ecology, single-species reaction-diffusion models generally deal with 
polynomial reaction terms. In a one-dimensional case, and if the environment is supposed to be 
homogeneous they take the form: 

t <") 

where u = u{t, x) is the population density at time t and space position x and D > is the diffusion 
coefficient. The function P, which stands for the growth of the population, is a polynomial of order 
> 1. The Fisher-Kolmogorov, Petrovsky, Piskunov (F-KPP) equation is the archetype of such 
models. In this model, we have P{u) = fiu — ^v? . The constant parameters ^ and 7 respectively 
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correspond to the intrinsic growth rate and intraspecific competition coefficients. In this model, 
the lower the population density u, the higher the per capita growth rate P{u)/u. More complex 
models can involve polynomial nonlinearities of higher order. Examples are those taking account 
of an Allee effect. This effect occurs when the per capita growth rate P{u)/u reaches its maximum 
at a strictly positive population density and is known in many species [HEKM]- A typical example 
of reaction term involving an Allee effect is |15l [T8\ [25] : 

P{u) = rn (1 — u) {u — p) , 

with r > and p € (0, 1). The parameter p corresponds in that case to the "Allee threshold" below 
which the growth rate becomes negative. 

In the previous examples, the reaction terms were assumed to be independent of the space 
variable. However, real world is far from begin homogeneous. In order to take the heterogeneities 
into account, models have been adapted and constant coefficients have been replaced by space or 
time dependant functions. In his pioneering work, Skellam j29] (and later, Shigesada, Kawasaki 
and Teramoto [28]) mentioned the following extension of the F-KPP model to heterogeneous envi- 
ronments: 

Here, the values of p{x) and 7(x) depend on the position x. For instance, regions of the space 
associated with high values of p{x) correspond to favorable regions, whereas those associated with 
low or negative values of p{x) correspond to unfavorable regions. As emphasized by recent works, 
the precise spatial arrangement of these regions plays a crucial role in this model, since it controls 
persistence and spreading of the population [U [H [ini ESI [24l [26l [27] . Models involving an Allee 
effect can be extended as well to heterogeneous environments, as in [13[ 125] . where the effects of 
spatial heterogeneities are discussed for models of the type: 

-Qi = + ^i^) ^ [(1 - ^)(^ - P{^)) + ^{^)] ■ (1-3) 

We also refer to [19J for an analysis of propagation phenomena related to a reaction-diffusion model 
with an Allee effect in infinite cylinders having undulating boundaries. 

In this paper, we focus on reaction-diffusion models with more general heterogeneous nonlin- 
earities: 

-Q^ = ^-Q^+Y^ pk{x)v!' +g{x,u), for t > 0, X € (a, 6), (1.4) 

k=l 

for some interval (a, b) in M. 

Since the behavior of such models depends on the precise form of the coefficients, their em- 
pirical use requires an accurate knowledge of the coefficients. Unfortunately, in applications, the 
coefficients cannot be directly measured since they generally result from intertwined effects of sev- 
eral factors. Thus, the coefficients are generally measured through the density u{t,x) [3^. Prom 
a theoretical viewpoint, if u{t, x) is measured at any time t > and at all points x in the con- 
sidered domain, all the coefficients in the model can generally be determined. However, in most 
cases, u(t, x) can only be measured in some - possibly small - subregions of the domain (a, b) [35]. 
For reaction-diffusion models as well as for many other types of models, the determination of the 
coefficients in the whole domain (a, b) bears on inference methods which consist in comparing the 
solution of the model with hypothetical values of coefficients jl^, with the measurements on the 
subregions [31j. The underlying assumption behind this inference process is that there is a one-to- 
one and onto relationship between the value of the solutions of the model over the subregions and 
the space of coefficients. This assumption is of course not true in general. 
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In this paper, we obtain uniqueness results for the coefficients Hk{x), k = 1,. . . ,N, based on 
locaUzed measurements of the solution u{t,x) of (jl.4p . The major differences with previous works 
dealing with comparable uniqueness results are (1) the size of the regions where u{t, x) has to be 
known in order to prove uniqueness, (2) the number of parameters we are able to determine, and (3) 
the general type of nonlinearity we deal with. 

Uniqueness of the parameters, given some values of the solution, corresponds to an inverse 
coefficient problem, which is generally dealt with - for such reaction-diffusion equations - using 
the method of Carleman estimates [5l |16] . This method provides Lipschitz stability, in addition to 
the uniqueness of the coefficients. However, this method requires, among other measurements, the 
knowledge of the density u{9, x) at some time 6 and for all x in the domain (a, b) (see p| ISj 171 [HI [36] ) . 
The uniqueness of the couple {u,fj.{x)) satisfying the equation (|1.2p given such measurements has 
been investigated in a previous work [8], in any space dimension. 

In a recent work, Roques and Cristofol ^23j have proved the uniqueness of the coefficient 
in ()1.2p when 7(2;) is known under the weaker assumption that the density u{t, xq) and its spatial 
du 

derivative — (t,xo) are known at a point xq in (a,b) for all t G (0,e) and that the initial den- 
ax 

sity u{0,x) is known over (a, 6). This result shows that the coefficient fj.{x) is uniquely determined 
in the whole domain (a, b) by the value of the solution u{t, x) and of its spatial derivative at a single 
point Xq. The present work extends this result to the case of several coefficients fJ,kix), k = 1, . . . ,N : 
given any point xq in (a, b) we establish a uniqueness result for the iV— uple (^i(x), . . . , fj.]y{x)) given 
measurements of the solutions u{t, x) of (II. 4p and of their first spatial derivatives in (0, e) x {xq}, 
starting with N nonintersecting initial conditions. 



2 Hypotheses and main result 

Let (a, b) be a bounded interval in M. We consider, for some T > 0, the problem: 



D-^ = ^^Xk{x)v!' + g{x,u), tG(0,r), a;G(a,6), 

k=l 

du 

aiu(t,a) — (3i — (t,a) = 0, t > 0, 
ox 

du 

a2u{t,b) + ^2-Q^{t,b) = 0, t>0, 



u{0,x) = u^{x), xG(a,6), 

for some G N*, and for - unknown - functions /i^ which belong to the following space A4: 

M. := {4^ £ C^''^{[a,b]) s. t. is piecewise analytic on (a, 6)}, (2-1) 

for some rj € (0,1]. The space C^'"^ corresponds to Holder continuous functions with exponent 77 
(see e.g. jl2|). A function ip S C^''^{[a,b]) is called piecewise analytic if their exist n > and an 
increasing sequence (Kj)i<j<„ such that ni = a, Kn = b, and 

n-l 

for all X G (a, 6), ij{x) = ^X[Kj,K,+i){x)ipj{x), 

i=i 

for some analytic functions (pj, defined on the intervals [kj, Kj+i], and where X[Kj,Kj+i) the 
characteristic functions of the intervals [Kj,Kj-^-i) for j = 1, . . . ,n — 1. In particular, if ^ G Ai, 
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then, for each x G [a, h) (resp. x G (a, 6]), there exists r = r^^ > such that ip is analytic on [x, x + r] 
(resp. [x — r,x\). 

The assumptions on the function g are: 

g{-,u) G C°'''([a,6]) for ah n G M, ^(x, •) G C^(M) for all x G [a, 6] and c/(-,0) = in [a, 6]. (2.2) 

We also assume that the diffusion coefficient D is positive and that the boundary coefficients 
satisfy: 

ai,a2,/3i,/32 > with ai + /3i > and Q2 + /32 > 0. (2.3) 

We furthermore make the following hypotheses on the initial condition: 

nO > in (a, 6) and u° G C2'^([a,6]), (2.4) 

that is is a function such that {u^)" is Holder continuous. In addition to that, we assume 
the following compatibility conditions: 

ai7xO(a)-/3i(nO)'(a) = and -D {u^ {a) = g{a,Q) if /3i = 0, 
a2uO(6)+^2(u°)'(6) = and -£> (uO)"(6) = 5(6, 0) if /32 = 0. ^'^ 

Under the assumptions (I2.ip - (l2.5p . for each sequence {fik)i<k<N S 7W^, there exists a time 
"^(Mk) ^ (0) +0*2] such that the problem ("P^^^.)) has a unique solution n G Ci'^/2 (['-'' "^Cfifc)) ^ ^]) 
(i.e. the derivatives up to order two in x and order one in t are Holder continuous). In the sequel, 
even if it means decreasing T^^^-^ in some cases and dropping the indices (fik) and u^, we only deal 

with values of t smaller than T so that the problem (^(^j.)) is well posed. Existence, uniqueness 
and regularity of the solution u are classical (see e.g. [21]). 

Our main result is a uniqueness result for the sequence of coefficients {fJ'k)i<k<N associated with 
observations of the solution and of its spatial derivative at a single point xq in [a,b]. Consider 

ordered initial conditions u- and, for each sequence (//fc)i<fc<Ar, let Uj be the solution of {V^^ 
Our result shows that for any e > the function 

{l^k)l<k<N ^ {ui{-,Xo),dUi/dx{-,XQ))i<i<N ' 

is one-to-one. In other words, we have the following theorem: 

Theorem 2.1. Let N G N*, ilJ-k)i<k<N O'f^d {fik)i<k<N be two families of coefficients in A4. 
Let {u^)i<i<N be N positive functions fulfilling \2.4\) and ^2. 5\) and such that u^{x) 7^ tij(x) for 

vP vP 

all X G (a, 6) and all i 7^ j. Let u-i and Ui he the solutions of ("^j^^^,)) c-'^^d ("^j^^^))) respectively, on 
[0,T) X [a, 6]. We assume that Ui and Ui satisfy at some xq G [a,h], and for some e G (0,T]; 

ni(t,xo) = nj(t,xo), 

du, , , dui, ^ for all t G (0, e) and alii & {!,■■■ , N}. (2.6) 
— (t,xo) = ^(i,xo), 

Then fik = p-k on [a, b] for a// G {1, • • • , N}, and consequently Ui = Ui in [0, T) x [a, b] for all i. 

The main result in |23j was a particular case of Theorem 12.11 A similar conclusion was indeed 
proved in the case N = 1 and for g{x,u) = —^u^. In such case, the determination of one coeffi- 
cient //i(x) only requires the knowledge of the initial condition and of (u(t, xq), 9M/9x(t, xq)) 
for t G (0,e). When N > 2, the above theorem requires more than the knowledge of the initial 
condition for the determination of the coefficients: we need a control on the initial condition, which 
enables to obtain measurements of the solution of (^(^j^,))) starting from different initial con- 
ditions. A natural question is whether the result of Theorem 12.11 remains true when the number of 
measurements is smaller than A^. In Section HJ we prove that the answer is negative in general. 





it^k))- 
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3 Proof of Theorem EOl 

For the sake of clarity, we begin with proving Theorem 12.11 in the particular case N = 2 (the proof 
in the case N = 1 would be similar to that of [23j, which was concerned with g{x, u) = —^v?'). We 

then deal with the general case of problems ('^(^j.)) (^(^^)) with > 1. 
3.1 Proof of Theorem EIH case = 2. 

Let and {fJ,2,fj-2) be two pairs of coefficients in 7W. Let ii5(x), ^2(2;) be two functions 

verifying (I2.4p and (12. 5|) and such that Ui{x) ^ U2{x) in (a, 6). Let ui and ui be respectively the 

solutions of {V^l^^^) ("P-^^ j^^) and U2 and U2 be the solutions of {V^l^^^) (^^1/^2)' 
We set, for z = 1, 2, 



Ui := Ui — Ui, mi := ^1 — jli and m2 := ^2 — /i2- 



The functions Ui satisfy: 



(3.1) 



for t € [0, T) and x G [a, 6], where 



bi{t,x) = fli{x) + /i2(a;) {ui{t,x) + Ui{t,x)) + Ci{t,x), 
g{x,Ui{t,x)) - g{x,Ui{t,x)) 



Ci{t,x) 



. du 

h{x, s) = s (mi(x) + m2(x)s) , 



Ui{t,x) - Ui{t,x) 
{x,Ui{t, x)) if Ui{t,x) =Ui{t,x 



if Ui{t,x) / Ui{t,x), 



(3.2) 



and the boundary and initial conditions: 



9^7- 

ai[/,(t,a)-/3i— ^(t,a) = 0, t > 0, 
ox 

dU- 

a2U,{t,b) + /32-^{t,b) = 0, t>0, 
[ Ui{0,x) = 0, xe{a,b). 



(3.3) 



Let us first assume that xq < b, and set: 

A+ = ^x > xo s.t. mi{y) = m2{y) = for all y G [j;oia;]|, 

and 

{sup(^+) if ^+ is not empty, 
Xq if A+ is empty. 

If xi = b, then mi(x) = m2{x) = on [xq, b]. Let us assume on the contrary that xi < b. 

Step 1: We show that there exist 6 > 0, X2 G {xi,b) and j G {1,2} such that the function 
{t,x) i-> h{x,Uj{t,x)) has a constant strict sign on [0,0] x (xi,X2], i-e. h(x,Uj{t,x)) > or 
h{x,Uj{t, x)) < in [0,0] x (xi,X2]. 
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To do so, let us define, for all x € [xi,b): 
mi{x) 



zix) 



m2{x) 



if m2(x) 7^ 0, 



^1 iy) (3 4) 

lim — — - if mi(x) = m2(x) = and m2(y) 7^0 in a right neighborhod of X, ^ ' ' 

y^x+ m2{y) 

+00 otherwise. 

mi{y) 



Whenever m2{x) ^ 0, z{x) is a root of the polynomial •). Notice also that the limit lim 



y^x+ m2{y) 

in the second case of the definition of z{x) is well defined since mi and ?Ti2 are analytic on [x, y] for 
y — X > small enough. 

Since, u1{xi) ^ i^{x\)^ we have — > for some j E {1, 2}. Moreover, there exists 

5 > such that xi + (5 < 6 and 

|ii^(x) — z(x)\ > r > on [xi,xi for some r > 0. (3-5) 

Indeed, if z{x\) = ±oo, we clearly have (|3.5p . If 2;(xi) / ±oo, from ()3.4p . z is continuous in a right 
neighborhood of x\. As the function -u*- is also continuous in this neighborhood, we get 



Moreover, Uj(t,x) E C^;^/2 (P'^) 

X [a, 6]) . This implies that Uj is continuous at (t, x) = (0, xi). 

As a consequence, there exists > small enough so that x) — u^(xi)| < r/4 in [0,^] x 

[xi, xi + 0], whence 

\uj(t,x) -u]{x)\ < ^ in (t,x) E [0,6^] x [xi,xi + 6l]. (3.6) 
Finally, setting 8' = min{9,5}, we have, from (j3.5p and (j3.6p : 

|uj(t,x) - z(x)l > I \uj(t,x) - n°(x)| - |n°(x) - z(x)| | > ^ > (3.7) 

for all {t,x) E [0,6*] x [xi,xi + 5']. 

Now, the definition of xi and the piecewise analyticity of mi and m2 imply that there exists 
6" E (0,5') such that: 

for all X E (xi,xi + 5"], the polynomial function /i(x, •) verifies h{x, •) ^ in M. (3.8) 

Indeed, assume on the contrary that there is a decreasing sequence y„ — )• xi such that h{yn, •) = 
in R. Since the functions h{yn,-) ^ire polynomial, we get mi(y„) = m2{yn) = for all n E N. 
Besides, as mi and m2 belong to for n large enough, mi and m2 are analytic on [xi,y„], which 
implies that mi = m2 = on [xi,x„] for n large enough. This contradicts the definition of xi and 
we get dHSD- 

From the expression (j3.2p of h and using (j3.8p . we observe that for each x E (xi, xi + 6"], either 
s = is the unique solution of h{x, s) = or m2(x) ^ and the equation h{x, s) = admits exactly 
two solutions s = and s = z{x). Let us set X2 = xi + 5" E (xi,6). From the strong parabolic 
maximum principle Uj{t,x) > in [0,6*] x (a, 6). Thus, using (13. 7p we finally get: 

h{x,Uj{t,x)) 7^ in [0,6*] x (xi,X2]. 

This concludes the proof of step 1. 

Step 2: We prove that xi = b. 
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From Step 1, let us assume in the sequel - without loss of generality - that (t, x) i-^ h{x, Uj{t, x)) 
is positive on [0,6] x {xi,X2] (the case h{x,Uj{t, x)) < could be treated similarly). Then, from 
the definition of xi, we deduce that h{x,Uj{t, x)) is nonnegative for all {t,x) S [0,0] x [xq,X2]. 

Since h{x2,Uj{0,X2)) > and Uj{0, •) = 0, it follows from ()3.ip that 



dt 



(0,X2) = h{x2,Uj{0,X2)) > 0. 



Thus, for e' £ {0,9) small enough, Uj{t,X2) > for t € (0,e'). As a consequence, and from the 
assumption of Theorem I2.H Uj satisfies: 



bj{t,x)Uj > 0, t G iO,e'), X € [xo,X2], 



dt dx"^ 

Uj{t,xo) = and Uj{t,X2) > 0, t G (0,e'); 
, Uj{0,x) = 0, X € {xo,X2). 

Moreover, the weak and strong parabolic maximum principles give that Uj{t,x) > in (0,e') x 
(xo,X2). Since Uj{t,xo) = 0, the Hopf's lemma also implies that 



dx 



{t,xo) > for all t E (0,e')- 



yi 



This contradicts the assumption (j2.6p of Theorem 12.11 Finally, we necessarily have xi = b and 
therefore mi = m2 = on [xo,6]. 

Step 3: We prove that mi = m2 = on [a, b]. 
Assuming that xq > a, and setting: 

A- = |x < s.t. mi{y) = m2{y) = for all y € [x,xo]|, 

and 

inf(^_) if A- is not empty, 

Xq if A- is empty, 

we can prove, by applying the same arguments as above, that yi = a and consequently mi = m2 = 
on [a, Xq]. 

Finally, mi = m2 = on [a, b] which concludes the proof of Theorem 12.11 in the case = 2. □ 

3.2 Proof of Theorem 12. IL general case > 1. 

We set for ah i,k e {I,-- - , A}, 

Ui := Ui - Ui, mk ■■= Hk - fj-k- 

The functions Ui satisfy: 

dUi ^d'^Ui 



D- 



bi{t,x)Ui + h{x,Ui{t,x)), 



dt dx"^ 
for t € [0, T) and x G [a, b], where 

if + g){x,Ui{t,x)) - {f + g){x,Ui{t,x)) 



(3.9) 



bi{t,x) 



d{f + g) 
du 



Ui{t,x) - Ui{t,x) 

x,Ui{t,x)) if Ui{t,x) = Ui{t,x), 



if Ui{t,x) / Ui{t,x), 



N 



h{x,s) = ^mfc(x)s'', 



k=l 



7 



xi :-- 



and 

N 

f{x,u) = ^pk{x)vl^. 

k=l 

Moreover, the functions Ui satisfy the fohowing boundary and initial conditions: 

dU- 

aiUi{t,a) - ^i-^{t,a) = 0, t > 0, 

a2U,{t,b) + /32^{t,b) = 0, t>0, ^^-^"^ 
[ Ui{0,x) = 0, xe{a,b). 

Let us set 

A-\- = |x > xo s.t. ruk = on [xq, x] for all € {1, • • • , iV}|, 

and 

sup(^+) if A+ is not empty, 
xq if A+ is empty. 

If xi = b, then for all k, mk{x) = on [xq, b]. Let us assume by contradiction that xi < b. As in the 
case = 2, we prove, as a first step, that there exist 9 G (0,T), X2 € ixi,b) and j € {1, . . . ,N} 
such that the function {t,x) i— >• h{x,Uj{t, x)) has a constant strict sign on [0,^] x {xi,X2], i.e. 
h{x,Uj{t,x)) > or h{x,Uj{t,x)) < in [0,0] x (xi,X2]. 

To do so, observe that, from the definitions of xi and M, there exists 6 > such that xi + 5 < b 
and all functions m-fc are analytic on [xi,xi + 6] and not all identically zero. Therefore, the integer 

p = max € N, mk{x) = o(^{x — xi)^ ^ as x ^ x|,^ for all A; € {1, ... , iV}| 

is well-defined. Furthermore, the function h can then be written as 

h{x, s) = {x - xiY X H{x, s) for ah (x, s) € [xi, xi + 5] x M, 

where 

H{x, s) = Mi(x) s + • • • + Mn{x) 

and the functions Mi, . . . , are analytic on [xi, xi + 6] and not all zero at the point xi (namely, 
there exists i G {1, . . . , A^} such that Mj(xi) ^ 0). In other words, the polynomial H{xi, •) is not 
identically zero. Since its degree is not larger than A^ and since H{xi,0) = and the real numbers 
Ui{xi), . . . , n^(xi) are all positive and pairwise different, there exists j £ {1, . . . , N} such that 

i7(xi,n0(xi)) /O. 

By continuity of H in [xi, xi + 5] x M and of Uj on [0, T) x [a, b], it follows that there exist 9 € (0, T) 
and X2 G (xi, b) such that 

H{x,Uj{t,x)) / for ah (t,x) G [0,9] x [xi,X2]. 

Consequently, 

h{x, Uj{t, x)) / for ah {t, x) G [0, 0] x (xi, X2]. 

The remaining part of the proof of Theorem 12.11 in the general case A" > 1 is then similar to 
Steps 2 and 3 of the proof in the particular case N = 2. Namely, we eventually get a contradiction 

dU • 

with the assumption that -^{t,xo) = for all t G (0, e), yielding xi = 6 and nik = on [xo,b] for 
all /c G {1, . . . , A^}. Similarly, = on [a, xq] for all A; G {1, ... , A^}. □ 
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4 Non-uniqueness results 



This section deals with non-uniqueness results for the coefficients (/Ufe) in (^(^j,)) under assumptions 
weaker than those of Theorem 12.11 These results emphasize the optimality of the assumptions of 
Theorem 12.11 

1-Numher of measurements smaller than number of unknown coefficients. 

We give a counter-example to the uniqueness result of Theorem 12. II in the case where the number 
of measurements is smaller than the number of unknown coefficients A^. 

Assume that the coefficients . . . , are constant, not all zero, and such that the polynomial 

N 

f{x,u) = f{u) = ^ ^kU^ 

k=l 

admits exactly — 1 positive and distinct roots ,zis[-i- Assume furthermore that ai = 

a2 = (Neumann boundary conditions) and that g = Then for each i = 1, ... — 1, zi is 
a (stationary) solution of ('P^j^^))- Consider a similar problem with the coefficients /ifc = t fj,/. for 
T 7^ 1 and k = 1, . . . ,N. Then, again, for each i = 1, N — 1, Zi is a solution of (^(^^^))- 
particular, assumption (j2.6p is fulfilled at any point xq G [a, 6] for k G {1, . . . , — 1}. However 
ifiiix), hn{x)) / {fLi{x), /iiv(x)). 

This shows that the determination of coefficients {pk)i<k<N requires in general A^ observa- 
tions of the solution of (^(^j.)); starting from A^ different initial conditions. 

2~Lack of measurement of the spatial derivatives. 



We show that if hypothesis ()2.6p in Theorem 12.11 is replaced with the weaker assumption: 

Ui{t,XQ) = Ui{t,XQ), for alH € (0, e) and all i E {1, • • • , A^}, (4-1) 

then the conclusion of the theorem is false in general. 

Let {lJ-k)i<k<N £ and assume that qi = Q2 = (Neumann boundary conditions). Let 
{ui)i<i<N satisfy the assumptions of Theorem 12.11 and assume furthermore that the functions 
and g{-,u) are symmetric with respect to x = (a -|- b)/2, i.e. 

/ (n-)l<j<Ar(x) = (n°)i<i<Ar(6- (x - a)) 11 ^ r .1 

< tor ail X fc [a, D\. 

[ 9{x,-) = 9{b -{x- a), •) 

Let /Ufc := /Ufc(6 — (x — a)) for all x E [a, b] and k £ {!,■ ■ ■ , A^}. 

Then, we claim that the solutions Ui and Uj of ("^j^^^)) ^-^d ("P^^^p satisfy (j4.ip at xq = and 
for e small enough. Indeed, we observe that, for each i G {1, • • • , A^}, Ui{t, 5 — (x — a)) is a solution 
of ('P(^j^))- By uniqueness, we have 

Ui{t, x) = Ui{t, b — {x — a)), for all t G (0, T), all x G [a, b] and alH G {1, • • • , A^}. 

In particular, Ui{t, ^) = Ui{t, ^) for t G (0,r) and « G {1, • • • , A^}. 

This shows that the assumption (j4.ip alone is not sufficient to determine the coefficients 

{fJ'k)l<k<N- 
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The above result is an adaptation of Proposition 2.3 in ^23j to the general case > 1. 
3-Time- dependent coefficients. 

We show here that the result of Theorem 12.11 is not true in general when the coefficients (/ijt) 
are allowed to depend on the variable t. 

We place ourselves in the simple case = 1 and gf = 0, and we assume that ai = 0:2 = 
(Neumann boundary conditions). We assume that (a, 6) = (0, vr), and we set u{t,x) = 1 + t cos^(x) 
and u{t,x) = 1 + t sin^ (23;) . 

Let us set 

1 / du 0'^u\ 1 / du d'^u\ 

i.e., 

, , UDt + 1) cos2(x) -2Dt , , , (16Dt + 1) sin2(2x) - 8Dt 

ai(t,x) = ^— ^ and Lii(t,x) = ^- — ^ . 

' l + tcos2(x) ' ' l + tsin^{2x) 

Then, for each t G [0, +00), /xi(t, •) and /ii(t, •) belong to AA. The functions u and u are solutions of 
('^(Vi)-' ^^■^'^ ("^(Ai))' respectively, and they satisfy the assumption (|2.6p of Theorem 12. II at xq = vr/2. 
However, the conclusion of Theorem 12.11 does not hold since fii ^ fli. 

^--Unknown initial data. 

We show here that the result of Theorem 12.11 is not true in general if the functions Ui and Ui 
are solutions of (^(^^)) a-^d ('P(^^))) with u^ ^ iq. This means that the coefficients (/x^) cannot be 
determined, given only the measurements {ui{t,xo),dui/dx{t,xo))i<i<N for t £ (0,e), if the initial 
conditions u^ are unknown. 

We build an explicit counter-example in the simple case = 1 and g = 0, with ai = 02 = 0. 
Assume that (a, 5) = (0, vr), and let us set u{t,x) = (1 + cos'^{x)) e^^ for some p > 0, and u{t,x) = 
(1 + sin2(2x)) e^^. We furthermore set 

1 / du d'^u\ 1 / Ou d'^u\ 

i.e., 

(AD + p)cos'^(x) + P-2D , (16D + p)sin^(2x)+p-8D 

'^^("^ = TT^) ^^^"^ = l + sin2(2x) ' 

thus pi and pi belong to Ai. Besides, u and u are solutions of (P^+™^'(")) and (pi+«-'(2.)^ 
respectively, and they satisfy the assumption (|2.6p of Theorem 12.11 at xq = it/2. However, we 
obviously have pi ^ pi. 

5 Numerical determination of several coefficients 

In the particular case A^ = 1, it was shown in |23j that the measurements ()2.6p of Theorem 12.11 
are sufficient to obtain a good numerical approximation of a coefficient pi{x). In this section, we 
check whether the measurements ()2.6p of Theorem 12.11 also allow for an accurate reconstruction of 
A^ coefficients {pk)i<k<N in the cases N = 2 and A^ = 3. 



10 



-6L 



0.2 0.4 0.6 0.8 

X 

(a) Determination of 2 coefficients 



0.2 0.4 0.6 0.8 1 

X 

(b) Determination of 3 coefficients 



Figure 1: (a) Plain lines: examples of functions fii (red line) and /i2 (blue line) in E. 
Dashed lines: the functions ii\ (in red) and ^2 (i^ blue) obtained by minimizing G(^^^^^^y In 
this case - /ii||i2(o,i) + ||/i2 - ^2IIl2(o,i) = 0-15, and Gi^^^^^^)[{jx[, ^i^)] = 9 • lO"^. (b) Plain 
lines: functions ^1 (red line), ^2 (blue line) and ^3 (black line) in E. Dashed lines: the 
functions /i| (in red), ^2 (i'^ blue), (in black) obtained by minimizing G(^^^^^^^^^y Here, 
11/^1 - MiIIl2(o,i) + 11/^2 -^2IIl2(o,i) + 11/^3 - /is ||l2(o,i) =0.38 and G(^^^^^^^^^)[{iil, ^i^, = 3-10-^. 



Given the initial data and the measurements Ui{t,xo) and -g^{t,xo), for t G (0,e) and 
i G {1, • • • ,N}, we can look for the sequence {fJ'k)i<k<N as a minimizer of some functional G{^^^y 
Indeed, for any sequence {jlk)i<k<N in A4^ , the distance between the measurements of the solutions 
Ui of (^(^j.)) and Ui of ^ ■ ) ^} can be evaluated through the function: 

G(^,,)[(/ifc)] = X] ll^i(-'2;o) -Ui{-,XQ)\\L2(Q^e) + ll"^(''^o) - -^(•,a;o)||L2(o,£). 

Then, G(^j,)[(^fc)] = and from Theorem 12.11 this is the unique global minimum of G(^^) in . 

In our numerical computations, we fixed (o, h) = (0, 1), D = 0.1, ai = a2 = and /3i = /32 = 1 
(Neumann boundary conditions), xq = 2/3 and e = 0.3. Besides, we assumed the coefficients ^k to 
belong to a finite-dimensional subspace E C A4 : 



E 



p e C0'''([0, 1]) I 3 {h,)o<j<n G M"+i, p(x) = ^hj-J ((n - 2) (x - c,-)) on [0, 1] \ , 

j=0 



(x — 2)^(x + 2)^ 

with Cj = ^^2 and J{x) = — if a; G ( — 2,2), and J{x) = otherwise. In our 

computations, the integer n was set to 10. 

Case N = 2 : 25 couples of functions {^1,112) have been randomly sampled in E"^ : for k = 1 
and k = 2, the components hj, in the expression 

n 

k^k{x) = Y,h]-J {{n - 2) (x - c^)) , 
i=o 

were randomly drawn from a uniform distribution in (—5, 5). 
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Starting from the initial data Ui = 0.1, and U2 = 0.2, the corresponding values of xq), 
^^{t,XQ), U2{t,XQ), ^^{t,XQ) were recordecn, which enabled us to compute Gj-^^^^j) [(/ii, /i2)] for 
any couple (/ii,/i2) in E"^. The minimizationci of the functions ^^ad to 25 couples (/ij,/i2), 

each one corresponding to a computed minimizer of the function G(^^^^^^y 

The average value of the quantity — Aii||L2(o^i) + \\fi2 — /"2IIl2(o.i)) over the 25 samples of 
couples {^i,fi2) is 0.25. The corresponding average value of G(^^^^2)(/^ii /^2) l-^ ' 10~^. Fig. [I](a) 
depicts an example of a couple (/ii,/X2) in E, together with the couple (/U*,/U2) which was obtained 
by minimizing G(^^^^^^). 

Case = 3 : in this case, the minimization of the function G(^i^^2,/i3) more time-consuming. 
We therefore focused on a unique example of a triple {fii, H2, Us) in E"^. The initial data were chosen 
as follows: Ui = 0.1, = 0.2, and Ug = 0.3. Fig. [I](b) depicts the triple (//i,^2,^3) in E, together 
with the triple (/i^, ^21/^3) obtained by minimizing G(^j^^2,/i3)- 

6 Discussion 

We have obtained a uniqueness result in the inverse problem of determining several non-constant 
coefficients of reaction-diffusion equations. With a reaction term containing an unknown polynomial 
part of the form J2k=i ^k{x) u^, our result provides a sufficient condition for the uniqueness of the 
determination of this nonlinear polynomial part. 

This sufficient condition, which is detailed in Theorem 12. 11 involves pointwise measurements of 
the solution u(t,XQ) and of its spatial derivative du/dx(t, xq) at a single point xq, during a time 
interval (0,e), and starting with nonintersecting initial conditions. 

The results of Section [4] show that most conditions of Theorem 12.11 are in fact necessary. In 
particular, the first counter-example of Section [5] shows that, for the result of Theorem 12. II to hold 
in general, the number of measurements of the couple {u,du/dx){t,xo) needs to be at least equal 
to the degree (A*") of the unknown polynomial term. 

From a practical point of view, such measurements can be obtained if one has a control on 
the initial condition. Nevertheless, since our result does not provide a stability inequality, the 
possibility to do a numerical reconstruction of the unknown coefficients /z^, on the basis of pointwise 
measurements, was uncertain. In Section [5l we have shown in the cases N = 2 and = 3 - which 
include the classical models (jl.2p and (|1.3p - that such measurements can indeed lead to good 
numerical approximations of the unknown coefficients, at least if they are assumed to belong to a 
known finite-dimensional space. 
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